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Abstract. It is known that real Non-Gaussian stable distributions are unimodal, not ad- 
ditive strongly unimodal, and multiplicative strongly unimodal in the symmetric case. By 
a theorem of Cuculescu-Theodorescu [5] , the only remaining relevant situation for the mul- 
tiplicative strong unimodality of stable laws is the one-sided. In this paper, we show that 
positive a— stable laws are multiplicative strongly unimodal iff a < 1/2. 



1. The MSU property and stable laws 

A real random variable X is said to be unimodal (or quasi- concave) if there exists a G M. 
such that the functions F[X < x] and F[X > x] are convex respectively in (— oo,a) and 
(a, +oo). If X is absolutely continuous, this means that its density increases on (— oo, a] and 
decreases on [a, +oo). The number a is called a mode of X, and might not be unique. A 
well-known example due to Chung shows that unimodality is not stable under convolution, 
and for this reason the notion of strong unimodality had been introduced in [8]: a real 
variable X is said to be strongly unimodal if the independent sum X + Y is unimodal for 
any unimodal variable Y (in particular X itself is unimodal, choosing Y degenerated at 
zero). In [8J Ibragimov also obtained the celebrated criterion that X is strongly unimodal 
iff it is absolutely continuous with a log-concave density. 

Proving unimodality or strong unimodality properties is simple for variables with given 
densities, but the problem might turn out complicated when these densities are not explicit. 
In this paper we will deal with real (strictly) stable variables, where very few closed formulae 
(given e.g. in [T4] p. 66) are known for the densities. A classical theorem of Yamazato shows 
that they are all unimodal with a unique mode - see Lemma 1 in [1] for the previously shown 
one-sided case and Theorem 53.1 in [T2] for the general result, but except in the Gaussian 
situation it is easy to see that stable laws are not strongly unimodal, because their heavy 
tails prevent the densities from being everywhere log-concave - see Remark 52.8 in [12J. 

Having infinitely divisible distributions, stable variables appear naturally in additive iden- 
tities, a framework where they hence may not preserve unimodality. Stable variables also 
occur in multiplicative factorizations as a by-product of the so-called M-scheme (or M- 
infinite divisibility), a feature which has been studied extensively by Zolotarev - see Chapter 
3 in [14] and also [TT] for the one-sided case. Another concrete example of a multiplicative 
identity involving a stable law is the following. Suppose that X, Y are two positive variables 
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whose Laplace transforms are connected through the identity E[e xx ] = E[e A " y ] for some 
a G (0, 1) (up to a reformulation one could also take a > 1). Then one has 



where Z a is an independent positive a— stable variable. In such identities, one can be 
interested in understanding whether the factorization through Z a modify or not some basic 
distributional properties. 

From the point of view of unimodality, it is therefore natural to ask whether stable vari- 
ables are multiplicative strongly unimodal, in other words, whether their independent prod- 
uct with any unimodal variable remains unimodal, or not. A quick answer can be given 
in the symmetric case because the mode of a symmetric stable variable is obviously zero: 
Khintchine's theorem entails then that its product with any independent variable will be 
unimodal with mode at zero - see Proposition 3.6 in [5] for details. However, in the non- 
symmetric case the mode of a stable variable is never zero (this is obvious in the drifted 
Cauchy case a = 1 and we refer to [H] p. 140 for the case a / 1), so that such a simple 
argument cannot be applied. The following criterion established in [5] Theorem 3.7, is a 
multiplicative counterpart to Ibragimov's theorem: 

Theorem (Cuculescu-Theodorescu). Let X be a unimodal random variable such that is 
not a mode of X . Then X is multiplicative strongly unimodal if and only if it is one-sided 
and absolutely continuous, with a density f x having the property that 

(1.1) t i — y /x( e< ) i s log-concave in R 

when X is positive (resp. t h-> fxi~ et ) ^ s log-concave in R when X is negative). 

With a slight abuse of notation, in the following we will say that a positive random variable 
is MSU if and only if (11.11) holds. Cases of multiplicative strongly unimodal, positive variables 
with mode at zero and such that (11.11) does not hold, are hence excluded in this definition. 
But these cases are particular, and by the above remark no more relevant to the content 
of this paper which deals with stable laws. Besides, a change of variable and Ibragimov's 
theorem entail the following useful characterization for positive variables: 

(1.2) X is MSU -<=>- \ogX is strongly unimodal. 

In particular the MSU property is stable by inversion and also, from Prekopa's theorem, by 
independent multiplication. Another important feature coming from (II. 2p is that the MSU 
property remains unchanged under rescaling and power transformations, which also comes 
from the obvious analytical fact that (11.11) holds iff 1 1— y Kie ait f x (K 2 e a2t ) is log-concave for 
some ai el and a 2l Ki,K2 G R*. Notice however that the MSU property is barely connected 
to the strong unimodality of X itself (several examples of this difference are given in [5]). 

In this paper we are interested in the MSU property for positive a— stable laws. For 
every a g]0, 1[, consider f a the positive a— stable density and Z a the corresponding random 
variable, normalized such that 



X 



d 



Z a x Y 



(1.3) 




A > 0. 
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Before studying the MSU property for Z Q , in view of (11.21) one must first ask if \ogZ a is 
simply unimodal. A positive answer for all a 6 ]0, 1[ had been given by Kanter - see Theorem 
4.1 in j9], who also deduced from [I] the decomposition 

\ogZ a = a^logbaiU) + (a — l)/a\ogL 

where L is a standard exponential variable, U an independent uniform variable over [0, tt], 
and b a (u) = (sin(cm)/ sin(w)) a (sin((l — a)u)/ sin(M)) 1_a . The random variable logL is easily 
seen to be strongly unimodal, but \ogb a (U) is not - at least for a = 1/2, see the Remark 
before Theorem 4.1 in [10]. This leaves the question of the MSU property for positive stable 
laws unanswered, and our result aims at filling this gap: 

Theorem. The variable Z a is MSU if and only if a < 1/2. 

To conclude this introduction, let us give two further reformulations of fll.ip in the positive 
stable case. The first one lies at the core of our proof, whereas the second one is probably 
nothing but a mere curiosity. Since f a is smooth, differentiating twice the logarithm entails 
that (11.11) is equivalent to the inequality 

(1.4) (x 2 ax) + xfjx))f a {x) < x\f a (x)f, x>0. 

If m a stands now for the mode of Z a and if x a = inf{x > m a , f'^{x) = 0}, then we know 
from (53.13) in [12J - an identity which is proved there for a certain class of positive self- 
decomposable distributions, but it can also be obtained for positive stable laws after taking 
the weak limit - and the discussion thereafter that (11.41) is true for any x G [0,x a ]. Hence, 
the MSU property amounts to the fact that it remains true for all x > x a . From (53.13) in 
[T2] we also know that (11.41) is equivalent to the positivity everywhere of the function 

x ^ I (f' a (x-y)f a (x) - f a (x-y)f' a (x))y- a dy, 
Jo 

but this criterion is not very tractable because of the long memory involved in the integral. 
Thanks to the Humbert-Pollard representation for f a which will be recalled at the beginning 
of Section 3, we finally mention that (11.41) is equivalent to 

/ „(-l) n X~( 1+an A (-l)n x -(l+*n)\ f (_l)" x -(l+an)\ 2 

> (1 + an) 2 ^— j — — > ^— j — : — ) < > (i + an )l_i , 

\i£ T(-na)n\ J T(-na)n\ J ~ } T(-na)n\ J 

a strange inequality which would plainly hold in the opposite direction if the terms of the 
series had constant signs. 

2. Some particular cases 

In this section we depict some situations where the MSU property can be proved or 
disproved directly, thanks to more or less explicit representations for the density f a or the 
variable Z a . First of all, in the case a = 1/2 it is readily seen from the known formula 
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that Z l/2 is MSU. When a = 1/3, Formula (2.8.31) in [H] yields 

where Ky 3 is the Macdonald function of order 1/3, so that (11. lj) amounts to show that 
t i — y Ki/^(e t ) is log-concave. From (jl.4p and since i£i/3 is a solution to the modified Bessel 
equation 

x 2 K';„ + xi^ /3 = (x 2 + 1/9)JT 1/S 
over R, this is equivalent to (x 2 + l/9)K 2 ^ 3 (x) < x 2 (K^ 3 (x)) 2 for every x > 0. Because 
i^ 3 (x) < for every x > 0, the latter is now an immediate consequence of a Turan-type 
inequality for modified Bessel functions recently established in [1] - see (2.2) therein, so that 
Z 1/3 is MSU. 

With a probabilistic and more concise argument, one can also show that Property (11.11) 
holds for every a = 1/p, p > 2. A classical representation originally due to E. J. Williams - 
see Section 2 in [13J - shows indeed that after rescaling Z±, p is an independent product of p 
Gamma variables: 

z i/l = f Y ^iv x r 2 / P x ••■ x >Vi)/ P 

where each Yk/ P is MSU since its density is r(k/p)~ 1 y k ' p ~ 1 e~ y l{ y> o}- Hence, since the MSU 
property is stable by inversion and independent multiplication, it follows that Z\/ p is MSU 
for every p > 2. 

As a matter of fact, a much stronger property is known to hold when a = 1/p for some 
p > 2. In theses cases, it had namely been noticed in [TO] thanks to the classical computation 
of the fractional moments 

(2.D = 

for every s < a, and the duplication formula for the Gamma function, that the kernel 
(x,y) i — y f a (e x ~ y ) is totally positive - see pp. 121-122 and 390 in [10] as well as pp. 11-12 
therein for the definition of total positivity. In particular, it is totally positive of order 2 
(TP2) which means precisely that x 1— >■ f a (e x ) is a log-concave function - see e.g. Theorem 
4.1.9 in [10] Chap. 4. However, when a is not the reciprocal of an integer, the function 

T(l-s) 

S !->■ 



r(l-s/a) 

is not an entire function of the type - see e.g. [10] p. 336 for a definition - and by Theorem 
7.3.2 in [10], this entails that (x,y) i-> f a {^ x ~ y ) is no more a totally positive kernel. Karlin 
raised then the question whether or not it should be totally positive of some finite order - 
see [TO] p. 390, a problem which seems as yet unadressed. 

Let us finally show that Z 2 /3 is not MSU, in other words that (x, y) h-> f2/z{ eX ~ v ) is not 
a TP 2 kernel. Formula (2.8.33) in [14] (with a slight normalizing correction therein) yields 
first the expression 
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where Wi/2,1/6 is a Whittaker function. Hence, from formulae (6.9.2) and (6.5.2) in [6] we 
see that (jl.ip is equivalent to the log-concavity of t \-> g(e l ) with g(x) = e^U^x) and the 
notation U\(x) = \I/(l/6, A/3, x) for all A > 1/2, where 

is a confluent hypergeometric function (c > a > 0). We readily see that g'(x) = —e~ x U 7 (x) 
and g"(x) = e~ x U w (x), so that by (11.41) the MSU property for X 2 /3 amounts to 

{x 2 U 10 {x) - xU 7 (x))U 4 (x) < x 2 U^{x), x>0. 

Using twice the contiguity relation (6.6.5) in [6] and some simple transformations, we then 
find the equivalence condition 

(2.2) {xU^{x)-U 1 {x)/&){U 7 {x)-U A {x)) > -5U 2 (x)/Q, x > 0. 

Notice that (12. 2p is true as soon as x > 1/6 thanks to the obvious inequalities U 7 (x) > 
U±{x) > U\(x). However, some easy computations yield the asymptotics 

7K J r(i/6) ' v ; r(i/6) v ; r(5/6) 

when x — > + , so that (12. 2p does not hold anymore. This last discussion about the case 
a = 2/3 may look tedious, all the more that the proof that the MSU property does not 
hold for any a > 1/2 is quite simple as we will soon see. But this must be read as a 
preparatory example for the sequel, since we will need inequalities such as (12. 2p for some 
confluent hypergeometric functions in order to show the MSU property when a < 1/2, a fact 
which is more involved. 

3. Proof of the Theorem 

We begin with the only if part, and we will give three different arguments. The first one 
relies on the aforementioned Humbert-Pollard representation for f a - see e.g. formula (14.31) 
in 



f a ( x ) = \ sin(7rcm)r(l + an)x- {1+an) = ' ' 



re ^— ' n\ ^— ' r(— na)n\ 

n>l n>\ y ' 



X 



-(1+cm) 



Because a < 1 this expansion may be differentiated term by term on (0, +00), yielding 

1 

n>l 

and 



n>l v ' 

for every x > (these three series representations explain the reformulation of ( 11.41) in 
terms of a reverse Cauchy-Schwarz inequality mentioned at the end of Section 1). Using the 
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expansions up to n = 2 and the concatenation formula T(z + 1) = zT(z) we obtain, after 
some simplifications, 

2 -(2+3q) 

* 2 (^) 2 (*) - (**K(x)+xf a (x))f a (x) = 2r "_^ )r( _ 2Q) + o(x-( 2+3 ")) 

in the neighbourhood of infinity, which entails that (11.4j) does not hold when a > 1/2, since 
the leading term is then negative. 

The second one hinges upon an expansion for the density g a of the random variable 
Y a = logZ Q , which had been obtained in [2] - see (3.5) and (6.4) therein - independently of 
the Humbert-Pollard formula: 

9a(t) = e^- My E t b j a? +1 (-lYR j (-e- at ) 
j>o 

for every t e M, where the coefficients bj and Rj{x) are defined through the entire series 



1 = J2 b ^ and eZ+xeZ = J2(^r~ 



r(i + z) 



Besides, setting 

P a {x) = ^bja^i-iyR^-x) 

for any x > 0, we know from (6.3) in [2] that the series converge absolutely, and with exactly 
the same argument one can show that it can be differentiated term-by-term. On the other 
hand, simple computations entail that g a is log-concave over M if and only if 

P a {xf + x 2 {P' a {x)) 2 > P a {x){xP' a {x)+x 2 P':{x)) 

for every x > 0. Letting x — > + yields 

P a {xf + x\P' a {x)f ~ xP a (0) 2 = x ' 



a 



and 



rfl-a 



P a (x)(xP^(x) + x 2 P:(x)) ~ xP a (0)P^(0) = x ( — - - ) (l ' ' ' " " 1 



r(l-a)/ V r(l-2a) # 

where the evaluations of P a (0) and Pq(0) rest upon the definition of 6j and the fact that 
Rj(0) = 1 and R'j(0) = 2 J — 1. Similarly as above, we see that the second asymptotic is larger 
than the first one when a > 1/2. 

For the third and simplest argument, we will invoke an identity in law connecting two 
independent copies Y a and Y a of the random variable logZ Q , which can be readily obtained 
in changing the variable in Exercise 4.21 (3) of [3]: 

(3.1) Y a - Y a = U a 

where U a is a real random variable with density 

sin na 



u a (x) 



ir(e ax + 2cosvra + e~ ax ) 
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We compute then the second derivative of x t- > log(e ax + 2 cos7ra + e which is 

4a 2 + a 2 cos 7ra cosh ax 
(e ax + 2cos7r« + e~ ax ) 2 

Hence, we see that U a has a log-concave density over R iff a < 1/2. By Prekopa's theorem 
this entails that logZ Q does not have a log-concave density over K when a > 1/2, which 
means that Z a is not MSU. 

Remarks, (a) This negative result shows that the kernel (x,y) i— > f a (e x ~ y ) is not TP 2 when 
a > 1/2, which contradicts the affirmation made in Lemma 1 (iv) of [7J that this kernel is 
strictly totally positive for every < a < 1 (actually the contradiction could already have 
been seen in reading [TO] p. 390 carefully). Notice that this latter affirmation seems crucial 
to obtain the so-called bell-shape property for all a-stable variables with index a < 1 - see 
p. 237 in [7J. However, since this question is quite different from the topic of the present 
paper, we plan to tackle the problem (if really any) in some further research. 

(b) Though somewhat more technical, the methods resting upon Humbert-Pollard's and 
Brockwell-Brown's expansions give some insight on the location where the inequality fll.4j) 
breaks down, an information which could not have been obtained by the third argument. I 
had also believed for a long time that these two expansions would give the if part, but this 
still eludes me because of the alternate signs. 

(c) It would be quite interesting to see if the third argument could not give the if part either. 
From the analytical viewpoint this would be the consequence of a positive answer to the 
following question. If X is a real random variable with density such that the independent 
difference X — X has a log-concave density, does X have a log-concave density as well? 
This assertion, a kind of reverse to Prekopa's theorem for which we found neither references 
nor counterexamples in the literature, goes somehow the opposite way to the central limit 
theorem (which leads to a log-concave density after enough convolutions on any probability 
distribution with finite variance). Franck Barthe wrote me that he would be surprised if 
it held true in full generality. In our positive a— stable framework, the fact that \ogZ a 
is unimodal (and log-concave at both infinities when a < 1/2) might add some crucial 
properties, but overall I could not find any clue in the direction of this statement. 

We now consider the if part, using yet another argument since the three above methods 
turned out fruitless. We first prove two lemmas which are of independent interest. The 
second one is a generalisation of Williams' result and could be formulated in several ways 
(we chose the one tailored to our purposes). 

Lemma 1. Let X be a Beta (a, 0) variable andY an independent Gamma (c) variable such 
that (3 < 1 and a + (3 > c. Then the product X x Y is MSU. 

Proof: When a + /3 = c, the result follows easily without further assumption on /3 because 



(3.2) 



X x Y = Gamma (a), 
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a fact which can be found e.g. in [3] Exercise 4.2. When a + (3 > c, we first compute the 
density function of X x Y: two changes of variable entail 

T(a)T(/3)T(c) J u 
r(a + /3)x c ~ 1 f [°° -xu p-i, +1) c-(«+p) d , 
r(a)T((3nc)\ Jo 1 j 

(notice that this recovers (13.21) when a + = c), so that we will be done as soon as the 
function t !->■ g Qi/ 3, c (e*) is log-concave, with the notation 



Using now exactly the same discussion made at the end of Section 2 for the case a = 2/3 
(with adaptated computations) entails this log-concavity property is equivalent to 

(xg a ,p,c( x ) + (a + (3 - c)g a> p iC -i(x))(g aj/ 3 jC+ i(x) - g a ,p, c {x)) > ((3 - l)(g a ^^i(x)) 2 

for every x > 0. But in the above, the right-hand side is negative because j3 < 1, whereas 
the left-hand side is positive from the obvious inequality g a ,/3,c+i(x) > g a ,p,c( x ), an d since by 
assumption a + (3 > c. 



Lemma 2. For all integers p, n > 2 such that n > 2p, one has the following representation 
as an independent product: 

a n n 

Z~J n = — Beta(2/n, l/p-2/n) Gamma(l/n) x Beta (4/n, 3/p-A/n) Gamma(3/n) x ••• 

x Beta (2(p - l)/n, (p -l)/p- 2(p - 1) /n) Gamma((2j9 - 3) /n) 
x Gamma ((2p — l)/n) x •■• x Gamma ((n — l)/n). 

Proof. We first evaluate the fractional moments of Z^T using (12. ip . the duplication formula 
for the Gamma function - see e.g. Formula (1.2.11) in [6], and some rearrangement involving 
the crucial assumption n > 2p : for every s > —1/n one obtains 



E 



T{ns + l) r(s + l 

X 



r(s + i) r(ps + i) 

r(a + l/n) ... r(s + (n - l)/n) T(l/p) . . . T((p - l)/p) 



pp J T(s + l/p) ... T(s + (p-l)/p) r(l/n) . . . r((n - l)/n) 
r(a + 2/7i)r(l/p) \ / r(s + l/n) ' 



n" ' 



pp J \T(s + l/p)T(2/n)J V r(l/n) 
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X 



r(g + 4/n)r(2/ P ) \ / r(a + 3/w) 

r(a + 2/p)r(4/n) J V r(3/n) 

r( a + 2(p - l)/n)r((p - l)/p) \ /r(s + (2p - 3)/n) 



r(s + (p - i)/p)r(2(p - i)/n) ; V r((2p - 3)/n 

T(s + (2p-l)/n) \ x w /T(a+(n-l)///) 



x 



r((2p-l)/n) / V T((n-l)/n) 

On the other hand, it is well-known and easy to see that the fractional moments of the 
Beta (a, /3) and Gamma (c) variables are given by 

E [(Beta (a, /?))*] = ffi ± and E [(Gamma (c)) s ] = Q£±£L 

1 (s + a + p)L (a) 1 (cj 

The claim follows now by identification of the Mellin transform. |— I 

Remark. As mentioned before, we see from this proof that analogous product representa- 
tions for positive a-stable laws with any a rational can be obtained accordingly. This might 
be useful to some other problems. Compare also with Theorem 2.8.4 in [J4] where trans- 
forms of the densities f p / q are given as solutions to differential equations of higher order, an 
analytical representation which seems less tractable than Williams-type representations. 

End of the proof: We need to show (II. 4p for any a < 1/2 and x > 0. Setting g a (x) = 
(x 2 fa(x) + xfa(x))f a (x) — x 2 (f^(x)) 2 , we see from the Humbert-Pollard decomposition and 
its differentiations that the application a i— > g a (x) is continuous on (0, 1) for every fixed 
x > 0. By a density argument, it is therefore sufficient to prove fll.4j) for any a = p/n with 
p, n integers greater than two such that n > 2p, and every x > 0. This amounts to the 
MSU property for Z~? n , which now comes easily from Lemmas 1 & 2, the MSU property for 
Gamma variables and the stability of the MSU property by independent multiplication. 
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